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SUMMARY 


This study has been devoted to the development and evalu- 
ation of a theoretical model of the proprotor on a cantilevered 
wing, operating in normal cruising flight. This theory expresses 
the wing and blade motion in coupled form, and can include any 
number of mode shapes required to describe the motion accurately. 

It has been applied to the investigation of the dynamic character- 
istics of the Bell and of the Boeing design. The Bell rotors are 
gimballed and the Boeing rotors are hingeless. The analysis in- 
cludes the frequency response to gusts and cyclic pitch, and an 
eigenvalue analysis of the dynamic system. 

Based on the theoretical results included in this study, the 
following conclusions may be stated: 

(a) The choice of mode shape (rigid-body mode or elastic-coupled 
mode) affects the damping significantly. The dependency of 
the damping on the mode shape can be estimated for the first 
beam bending mode. The blade inplane deflection opposing 
the rotor direction of rotation, accompanied by the forward 
out-of-plane deflection, increases the damping, comparing it 
with the rigid-body calculation. The inplane deflection pro- 
ceeding in the rotor direction of rotation decreases the 
damping. The mode shape has little influence on the frequencies 
of the system. 

(b) The results of the frequency response are quite similar to those 
of Johnson, in spite of the difference in the mode shapes. The 
amplitude of the response is slightly different, since struc- 
tural damping was not included in the present calculation, and 
the mode shapes used were different. 

(c) The analysis of the eighteen degree-of- freedom system showed 
that the higher- frequency degrees of freedom have small influ- 
ence on the basic degrees of freedom. 
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Stresses or bending moments of the wing or blade can be 
predicted from the motions of the wing and blade obtained from 
this analysis. In addition, this analysis may be applied to the 
development of an automatic control device to alleviate the gust 
response of the vehicle. 
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SECTION 1 


INTRODUCTION 


1.1 General 

It has been recognized that the tilting proprotor aircraft, 
one of the composite aircraft family, is a very promising concept 
that combines into one aircraft the hover efficiency of the heli- 
copter and the high-speed efficiency of the fixed-wing aircraft 
(Refs. 1-3) . 

The typical tilting proprotor aircraft is a twin-engine air- 
craft with tilting rotors mounted on each wing tip. Its configura- 
tion consists of a fuselage, a high swept-forward wing, and an em- 
pennage. The empennage has a vertical stabilizer and rudder, and 
a horizontal stabilizer and elevator. The large diameter rotors 
are three bladed. hincreless or gimbal-type rotors which are mounted 
on the rotor shaft. The rotor shaft is connected through the gear- 
box to each engine in the pylon attached at the wing tip. The con- 
version system provides the rotation of the rotor pylon from the 
vertical position to the horizontal position and return, in order 
to obtain the helicopter mode or airplane mode corresponding to the 
desired flight regime. 

When the aircraft takes off or lands, the rotor pylon is ro- 
tated to the vertical position to achieve vertical takeoff or land- 
ing similar to the helicopter. The flight controls apply pitch 
changes to the rotor to provide the longitudinal and directional 
control corresponding to helicpoter rotor, cyclic pitch, while the 
collective pitch controls vertical flight and roll motion. 

In high-speed flight, the rotor pylon is rotated to a hori- 
zontal position similar to that of the conventional propeller type 
aircraft. The thrust is produced by the rotor, and the lift by the 
wing. The flight controls are provided by the conventional aircraft 
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control surfaces such as the elevator, rudder and aileron. 

In addition to the above modes, various conversion inodes 
can be obtained. At that time the rotor pylon tilts to some posi- 
tion between the vertical and horizontal, where it can be safely 
locked. This makes STOL-type operations possible. 

The tilting proprotor is exposed to a severe aerodynamic 
environment including gusts, the wake of preceding blades, and 
harmonic air loading like a helicopter. But its dynamic and aero- 
elastic characteristics are in many ways unique; for example, the 
large flexible blades with a large amount of twist experience sig- 
nificant coupled out-of-plane (flapping) and inplane (lagging) 
motion. 

As described later in Subsection 1.2, several years of ex- 
perimental and theoretical analyses have been conducted to estab- 
lish a fundamental understanding of the dynamic and aeroelastic 
behavior. However, it is necessary to understand the aeroelastic 
response of this aircraft to atmospheric turbulence more adequately 
and to predict it more accurately, since during the preliminary 
design phase, vibration level prediction is required in order: (a) 
to evaluate the fatigue life of the blade and wing, (b) to estimate 
the ride qualities of the vehicle, and, if necessary, (c) to develop 
suitable gust alleviation devices. 

Several design compromise concepts, which make the present 
analysis distinct from helicopter aeroelastic analysis, are now 
stated briefly. 

In order to obtain high hover efficiency from the rotor, it 
is desirable to achieve low disc loading, in other words to use 
large-diameter rotors whose swept discs reach nearly to the fuse- 
lage. When the aircraft is operated in high forward speed axial 
flight in the airplane mode, the rotor is operating at a high inflow 
ratio (the ratio of axial velocity to blade tip speed) . This pheno- 
menon is very different from the helicopter rotor operation which 
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involves low inflow. High inflow operation requires a large built- 
in angle of twist for efficient cruising. Therefore, significant 
coupled out-of-plane (flapping) and inplane (lagging) motion occurs 
in the large, flexible and twisted blade. This phenomenon makes 
analysis more complicated. 

The engines and gearboxes are usually located at the wing 
tip to avoid transmitting high power through a long drive shaft. 
This leads to low wing natural frequencies and possible resonances 
in the low frequency range. Also, the center of gravity of the 
pylon and rotor does not usually coincide with the elastic axis of 
the wing. Hence, this results in coupled bending and torsion. 

1.2 Brief Survey of Past Work 

A good review and an elementary description of the dynamic 
and aeroelastic problems associated with the tilting proprotor air- 
craft were given by Reed in Ref. 4 and Loewy in Ref. 5. In this 
report no attempt will be made to repeat the reviews given in Keis. 
4 and 5. The only references cited here will be those pertinent 
to the problem being treated. 

The possibility of propeller whirl flutter — a dynamic in- 
stability that can occur in a flexibly mounted aircraft engine- 
propeller combination — was first recognized in the late 1930's 
by Taylor and Browne (Ref. 6). It was not until 1960 that it be- 
came a problem of practical concern -- with the appearance of the 
turboprop aircraft. 

Following the two fatal turboprop aircraft accidents, it 
was established that propeller whirl flutter could have occurred 
if the nacelle stiffness was severely reduced, for example by a 
structural failure. Several generalized studies were conducted 
at NASA-Langley . One of them was carried out by Houbolt and Reed 
in Ref. 7; it gives an elementary and basic treatment of the equa- 
tions of motion and propeller aerodynamics for propeller-nacelle 
whirl flutter. 
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Because VTOL configurations have unconventional propeller- 
rotor systems, whirl flutter was a major design consideration on 
present proprotor aircraft. 

The analysis presented in Ref. 8 is for a two-bladed rotor 
free to tilt on a shaft with two nacelle degrees of freedom (pitch 
and yaw) . No lag or coning degrees of freedom are considered. 

The analytical method was compared with test results for an exist- 
ing tilting proprotor aircraft (the Bell XV-3) and of subsequently- 
tested scale models. They showed good agreement. 

Young and Lytwyn in Ref. 9 present a very precise analysis 
for the whirl stability of a multi-bladed rotor mounted on a nacelle 
which has pitch and yaw degrees of freedom. Each blade has one flap- 
wise degree of freedom. The blade mode shape is assumed to be a 
rigid body mode shape. It was concluded that whirl stability is 
poorest when the nacelle pitch frequency equals the nacelle yaw 
frequency, but in this situation nacelle damping is quite effective. 
There is an optimum value of flap bending frequency somewhere between 
1.1 auJ for highly sLabilized whirl motion. 

This analysis neglects (as do Refs. 4 and 8) the effect of 
coning on proprotor aerodynamics, and flap bending mode shapes other 
than the rigid blade mode used. Also, autorotation flight must be 
considered as well as powered flight. 

In Ref. 10, Gaffey points out. that a highly coupled blade mode 
has substantial flap bending even if the primary mode involves in- 
plane motion. This occurs in the case of a highly twisted blade or a 
blade operating at high geometric pitch angles such as a proprotor 
blade. The analysis shows that a moderate amount of negative 6^ 
(flapping angle at the blade root gives the pitch angle reduction of 
the amount 6 'tan if $3 is positive) has a stabilizing influence 
on proprotors subject to flap-lag instability at high inflows. 

Preliminary design studies of prototype vehicles (Refs. 11 and 
12) as a part of the current NASA/ARMY sponsored tilting proprotor 
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research aircraft program give some results from dynamic and aero- 
elastic analyses done by Bell and Vertol. 

Johnson (in Refs. 13 and 14) derived the equations of motion 
for a cantilever wing with the rotor at the wing tip. He develops 
a nine degree-of-freedora model which involves blade flapping motion 
and lagging motion (each has one collective and two cyclic motions, 
respectively), wing vertical bending, chordwise bending, and torsion. 
This model is applied to two proprotor designs and compared with the 
results of some full-scale wind tunnel tests. It shows reasonable 
correlation between theory and experiment. 

In conclusion, it appears that most of the investigative work 
has been concerned with whirl flutter. The above review shows that 
more knowledge is needed for the solution of tilting proprotor air- 
craft dynamic and aeroelastic problems. 

1,3 Objectives of the Present Study 

The objective of this study is to establish a verified method 
of predicting the dynamic and aeroelastic behavior of the tilting 
proprotor aircraft in order to evaluate the fatigue life of the 
blades and wings and also to estimate the ride quality. 

The equations of motion for a cantilever wing with a rotating 
rotor at the wing tip will be derived as consistently as possible. 

The great complexity of rotor blade motion will be included by ac- 
counting for blade rotation (i.e., centrifugal and Coriolis forces), 
significant inplane motion, and the large twist and high pitch angles 
at high inflows. 

The resulting system of equations, obtained using modal 
analysis, will be applied to the analysis of the two proprotor de- 
signs (one is a rigid, soft-inplane type rotor and the other is a 
gimballed, stif f - inplane rotor) . 

Finally, the eigenvalues and frequency response of each pro- 
protor design will be determined to establish their dynamic charac- 
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teristics; the results of the analyses will be compared in terms of 
eigenvalues for normal modes (coupled elastic inodes) , assumed modes 
(uncoupled elastic modes), and rigid body modes, using Galerkin's 
method. 
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SECTION 2 


THE EQUATIONS OF MOTION 

2.1 Model and Coordinate System for the Analysis 

The primary interest of this study is in the dynamic and 
aeroelastic phenomena of the wing, pylon, and rotor of the tilt- 
ing proprotor aircraft in cruising flight. Hence, the dynamical 
system considered here consists of a cantilever right wing with 
a pylon at the wing tip, and an N-bladed rotor mounted on the 
pylon, as shown in Fig. 1. The model will be restricted to the 
cantilever wing, since it is sufficient to obtain a basic under- 
standing of the proprotor motion, and many such proprotor models 
have been tested in wind tunnels. 

Therefore, the aircraft rigid body motions are neglected 
and the wing anti symmetrical modes are also dropped. The left 
wing motions including the pylon and rotor are given by the 
mirror image of the right wing. 

The wing is assumed to have a high aspect ratio, so that 
strip theory is used for the wing aerodynamics and beam theory 
for elastic bending and torsion. Wing sweep and dihedral will 
not be considered, but angle of attack and angle of twist {built- 
in twist) will be considered. The elastic axis is assumed to be 
a straight line. The elastic axis coincides with the y-axis as 
shown in Fig. 2. The free stream vector coincides with the z- 
axis. Therefore, the angle between the z-axis and the chordline 
of the wing results in a wing angle of attack (positive nose up) . 
The wing motion (Fig. 2) consists of elastic bending and elastic 
torsion. The deflection u w of the wing elastic elastic axis per- 
pendicular to the y-z plane is called vertical or beamwise bend- 
ing (positive upward) . The deflection w w parallel to y-z plane 
is termed chordwise bending (positive forward) . Torsion p w is 
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defined as pitch angle change (positive nose-up) . 

A pylon of large mass and moment of inertia is assumed to 
be rigidly attached to the wing tip. Therefore, the pylon motions 
(Fig. 2) corresponding to the wing motions are defined as vertical 
displacement r x along the x-axis (positive upward) , longitudinal 
displacement r z along the z-axis (positive forward) , pylon yaw v 
about the x-axis (positive counterclockwise) , pylon pitch v p about 
y-axis (positive nose-up) , and pylon roll v r about the z-axis (posi- 
tive for clockwise rotation, when looking forward) . The above py- 
lon motions are accounted for at the point where the wing elastic 
axis crosses the plane, parallel to x-Z plane, which includes the 
rotor shaft. The pylon lateral displacement along the y-axis is 
neglected as a higher order effect. 

The pylon and the rotor shaft are assumed to be parallel to 

the free stream in equilibrium flight (r , r , v , v , and v are 

x z y p it 

all zero) , regardless of the wing angle of attack. 

The rotor is located at the distance h fFig. 7 ) from the 
wing tip elastic axis to the rotor-hub (positive forward from the 
elastic axis) . The distance h is termed the mast height. The rotor 
consists of N blades, whose rotational speed SI is defined as positive 
if clockwise looking foward. The blade (Fig. 2) has out-of -plane 
(flapping) deflection w r , defined positive for forward displacement 
from the disc plane (upward in helicopter mode, while the rotor 
shaft is vertical), and inplane (lagging) deflection v , defined 
positive for clockwise deflection regardless of rotor direction of 
rotation. The lower case letter n means the nth blade, n=l,2,...N. 
Blade torsion is neglected here. 

The azimuth position T (Fig. 3) of the nth blade is de- 
fined as: 
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( 2 . 1 ) 


= at + + a 4 * r 

where f- is measured from the vertical and t is time. The 
n 

phase angle between blades, Aip n is defined as 
Afk* 5; CH -0 - 3 JL 

N 

The azimuth position of the nth blade excluding the pylon roll 
motion is denoted as 

= At + (2 * 3) 

2.2 Basic Formulation for Powered Flight 

The governing linear equations of motion are derived in 
this subsection. A more complete and detailed derivation is 
given in Appendix A. 

The equations involve ten unknowns: v n , w n » r^, r z ' v y f v p' 
v ,u , vr , and p {actually v and w represent N unknowns, res- 
pectively, but for convenience they may be treated as one unknown 
without inconsistency) . The equations consist of three categories 
the blade equations, the wing equations, and matching conditions 
between rotor and pylon. 

The blade equations are complicated by the pylon motions 
which produce the centrifugal forces and Coriolis forces. The 
wing equations are derived from beam theory. The rotor is rigid- 
ly attached to the pylon for the powered flight case and, there- 
fore, rotor motion is related to pylon motion; this gives the 
matching conditions between the rotor and the pylon. 

(a) Blade Equations 

-^[{ceD^o. + Ce« b c~’e.}i£L] 

+ -(«)„) ^ a* '~ ,0 » 
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-■§r( T '! ! r") + m{ w. +r, + 

- r 1* o«> 

+ rij Ain $*] = (P^ (2.4) 


cod 2 ©a +C^X)e^ a ^. a ©» 


■^3 [{< E »c 

+ "fr* ***■ ®i> '^^■j 


1 aV w ~ 

J sr 1 

* 


3r 


( T + ™ { V, - jO? v* - r, Ain <f „ 


+ - A i>p "1 <£„ - 6 py 4^° §H } 

= (Pain 


( 2 . 5 ) 


wnere 

r 


running spanwise coordinate for the 
blade from the axis of rotation. 


6 angle formed by the rotor disc plane 

and the blade sectional chordline, 
usually including built-in angle of 
twist and collective pitch. 

(EI) c bending stiffness in the blade sec- 

tional chordline direction. 


(EI) b 


m 


bending stiffness in the direction 
perpendicular to the chordline. 

spanwise mass of the blade per unit 
length . 
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T 


P 

z 


p 


e 


( ) 


n 


centrifugal force at r expressed as 

rft 




m r dr 


resultant force per unit length in 
the z direction on the blade (posi- 
tive forward) . 


resultant force per unit length in 
the circumferential direction on 
the blades (positive clockwise when 
looking f owardT . 


about nth blade. 


The blade geometry is shown in Fig. 3. 


(b) Wing Equations 


-|lj jj(ei„) c 4^J 

+ -|^[{(Eiw) e -(«.),] ^ a. 

+ m w U w +S«fv = F* u.6) 

+ cei-w)® 



( 2 . 8 ) 


3 

*2 


■{cw-fM + iwi + s. = m 3 


where 

e w 


(El ) _ r (El ) r. 

W C W b 


(a V 

m w 


w 


a 


angle between the wing sectional chordline 

and free-stream direction. 

wing bending stiffness defined similarly 

as for the blade. 

wing torsional rigidity. 

spanwise wing mass per unit length. 

wing mass moment of inertia about the 

elastic axis per unit length. 

static mass moment of wing segment about the 

elastic axis defined as 

s - ms 
a 

where s is a distance between the center of 
gravity and the elastic axis of the wing 
and positive if center of gravity is ahead 
of the elastic axis. 

F , F_ resultant wing force per unit length in the 

X z 

x and z directions , respectively. 

My resultant wing pitching moment per unit 

length. 

The wing cross-sectional geometry is shown in Fig. 4, and 
pylon motions are expressed with wing deflections as 




' Uw 1 

Vz 




‘ =, 

st y 

u ? 


1 fw 

vu 



^ * r ) 

1 

J 


(2.9) 


9*t 
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(c) Matching Conditions between the Blade and Pylon 

Mi r x - -|y [|cei») c 6 w + Ce i*) b j 

j^(EX*) c “ (EXw)^| ,o<**.'Qa Cv*> Ov 

+ £ [J*w{>i + k'C> f -C-m, *i«$„)jar] 

-t [(Va^fc.+*B 0<4 * (2] 

H*| t Jo J* 


M f »i - 




.[{cei w )« «™*e«r + CeiJ, ^4.} 


2 

^3 


j{(EI,) e -(EI,)»J ^dm^Sw -^|rj 
Ilf ( + vtT«) 4 r] 




=Z[j*(H.B^€H + ^ft ce »?K +Pr)arj 

1*2 * + [{^lw) c K E U, A^dwj j 

+ [{cej,^ - ( ei w .),j ^ a, c«, e, -§^? j 

+r[f„ w »{(* 1 + T> «J + nr“^ + r(w„Zi«<EJ 

- 2 £X r ( «;« l#> i J - <t ( v* aw $ jJ d r j ^ 

= - l i&) + + rPz./^K$, 


(2.11) 
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+ Pi V, cco®, - uT„ Pe coo j dr 

Ip, ^ + [U) -ag L +Z [\\{- *( u„ ** o) 

+ 4 (r x +<^ P ) - 2 ar -rCur w coo$ H ) 
-n^] dr]^ = Z 
+ P z ) - vj;p e •+■ \r H P z >**-<- ^* 

-irP z coo$ M Jdrj 


( 2 . 12 ) 


(2.13) 


Vi - [I cei ^c +^6w + (EIw)b ^ aQ ^} “I 

~ [{CEI W ) c ~ C=Iwib} ^ Rr 000 'J 


iK 


'M(r*LL + r tr H ) ar 


+ S,fC 

■fi[C‘'«*L 


1 


where 


M Pylon mass 
P 


I Pylon yawing mass moment of inertia 

p y 

I Pylon pitching mass moment of inertia 

P P 

I Pylon rolling mass moment of inertia 


(2.14) 


(”),{*) Time derivative applied to entire formula 
in the parentheses 
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The above equations are derived independently of Refs. 7, 

9, and 13. However, if the wing motions are eliminated, the 
equations are similar to the basic whirl flutter equations derived 
in Ref. 7 or Ref. 9. 


Supplement for Autorotational Flight 


Autorotation may be defined as the condition of flight where 
there is no restraint of the rotor rotation about the rotor shaft. 
Therefore, no rotor torque is transmitted to the shaft, and no pylon 
roll motion is transmitted to the rotor. This means that rigid body 
rotation of the entire rotor about the shaft will be produced. This 
rigid body rotation is designated as v R (positive clockwise). 


The equations of motion for autorotational flight are almost 
the same as those for powered flight. The equations: 2.4, 2.6, 2.7, 
2.8, 2.10, 2.11, 2.12 and 2.13 are the same. In Eq. 2.5 V r must 
be replaced by v_. In lieu of Eq. 2.14 the fact that the rotor 
inplane motion is independent of the pylon roll motion in auto- 
rotation flight results in Eq. 2.15 as follows: 


EflW 1 * -V rU-nWH 

•Mel l '<* -U 

= C( JVp.) ar] 

This yields a new degree of freedom. 

The angle in Eq. 2.1 is also expressed as 

- nt + vu + 

in this case. 
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SECTION 3 


AERODYNAMIC FORCES 


3.1 Rotor Aerodynamic Forces 


The rotor aerodynamic forces will be evaluated next. The 
analysis is almost the same as in the helicopter blade element 
theory. The significant difference is that it is impossible to 
assume small angles for the angle of twist, the collective pitch 
and the angle of attack of the blade operating in the high inflow. 
The basic idea is presented by Young and Lytwyn {Ref. 9). 

The section aerodynamic lift and drag forces dL, dD yield 
the resultant forces P and P Q as follows (in Fig. 5} , 

Z t7 


P z = dL cos - dD sin <f>^ 
P* = -dD cos <J>. - dL sin <f>. 

U 1 1 

where 

dL = 1/2 pac B U 2 a 
dD = 1/2 p CgD 2 ^ 


(3.1) 


(3.2) 


and the angle <J>^ is defined as. 


sin *i 

|u I (3 * 3 

cos*. = I U tI 

u 

The blade section inplane velocity (positive counterclock- 
wise direction) is and Up is the blade section out-of-plane 
velocity (positive for negative z direction) . Therefore, the re- 
sultant air velocity is expressed as, 


U * J uj + 

In cruise flight, 

U T - ttr 

U = V + v 
P 


(3.4) 

(3.5) 


16 



where V is the forward aircraft speed (.axial velocity) and v is 
the induced inflow velocity. The induced velocity is very small 
compared with the forward velocity V in the high inflow operation 
Therefore the induced velocity is neglected in this entire study. 

The inflow ratio X is defined as the ratio of the axial 
velocity to tip speed of the blade: 


V 
UR 

The section effective angle of attack is 


(3.6) 


4 Opp 


— "d x*0 l 



+ 0 C 


(3.7) 


the angle 0 g formed by the disc plane and the blade chordline in- 
cludes the angle of twist and collective pitch. Therefore 


where 0 ftT is the built-in section angle of twist of the blade, 
which is zero at 75% span of the blade, and 0 ?5 is the collec- 
tive pitch to be obtained from the performance calculation. 
Usually 0 is: 


0.75 



-V ) 

/ 


4 * 


(3.9) 


The first term express the inflow at 75% span and depends on 
the proprotor design to obtain the optimum cruising performance. 
Pitch-flap coupling and blade pitch control 9 will be dis- 

cussed later. Finally: 


tf 


4“ 0.7S + Qpf 



+ Qt 


(3.10 
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Substituting Eqs. 3.2 and 3.3 into Eq. 3.1 yields 


2 X = -i.fc 5 a U|Ur|c* --jPCbCdoUU* 

P»= fi.[~i : /-c e C..U|U T | --Lf>%a UUt«] (3-U) 

where 

sign (Cl) = -5^. 

I-TlI 

= Si (3.12) 

Next the perturbation method is applied to Eq. 3.11 to derive the 
aerodynamic forces of the rotor in disturbed motion: 

& = CP,). + (-t%)> T 4 (-§&). SUp 4 & 

P* = CP.)o + + (-*£).&* + (#)> <3 - 13 ’ 

In these equations, ( ) means those values are evaluated 
in the equilibrium state, given by the following expressions: 

(&)<,= iP c * a cnRj*[«.T, -%sXr» 

-4- A a oc<> T, - A 5 Toj 

(H£l - -k fC * tt WR [ 2 ** Tl + * (| - %-)•*. fi 

(■Hi=7<’ c * a|n|R [-0 + %-) r ‘ + *** - 2 ^H 

1 L (3.14a) 
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( a). - t f =» a (aR? 

v, - xV t.] n. 

(■f^ a -T' ,c » a|a|R [- 2 -% iT »- Xo< ' T ' 

^-|BlJ =i^c B aln|R[-a.T 2 + x(l--^-)r, 

(-flrl = l f c,a(QR) ! [- A t* - A J *] SI (3 - 


where 


9C.^ 

T « = 7 = = =~ (4 = 0 ; 1 , 2 , 3 ) <3.15) 

J * + A 


oc = 

r A 


(3.16) 

X = 

1 A 1 


(3.17) 

<Xo = 

Oat +0-TS 

+ -i Up 

" W 1 UA 

(3.18) 


The perturbation quantities 6 u t , and 6a will be con- 

sidered next. The perturbation velocities 5U T and 6Up consist of 
rotor and pylon motion and gust velocities. 

The velocity at the position of the blade axis in x,y,z 
coordinate (Fig. 2, see details in Appendix A) can be written as 
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^ - n ^ AM*, - £1 'U'n ^ ***< + ^X 

“Ht^p - V r - O.r0r coc 

5 - Xlrcoo^K + rtV coo</^ -0.rU r 

«»</•*- _C1 V K ^K. </* 

Z= VJ^ +?k - r b> p coo <k 

+ ri^^M.4 K +Xlr^coo^ 


(3.19) 


Wind velocities at the blade element in the disc plane U T are then 
obtained by applying the transformation matrix due to pylon motion 
to change the x,y,z coordinate system to blade disc coordinate 
system. Wind velocities 'U xB and U yfi are x B , y g components of the 
inplane velocity of the blade in the x B ~y B coordinate system fixed 
to the rotor hub (Fig. 6). Gust velocity as well as the aircraft 
forward velocity is included. 


f u»j 


] j 

* — - 

■ 

l U P J 



I 

-*V 


-V’r. 


-u, 



Tangential velocity U T is expressed as 

= ~ U xa 4*. + 

After higher order terms are neglected, 

U T = ril - ^ coo</ K 


(3.20) 


(3.21) 






C<jO</ K ) 


20 



(3.22) 


u P - v-^ 

+ ^+ h + r 4- k 

Finally the perturbation velocities SU^ and 60^ are given as 

Su T » Ixl|r[-a^^<*, +X-^c 

4 .) + A (if. + «.<£,)] 

Sut*|n|R[^-^ +-^* + -£■ -x\&4, 


o 

_ 

R 


+ lJ, c^o ( 

7 


where 


V ’ R * R 

+ i^A<f K ] (3.23) 

V J “ SCI A|1>) 

(3.24) 

Pitch-flap coupling gives a change in blade pitch angle propor- 
tional to the flapping angle at the root. It is defined as 


(a ) 3*% 

^k—^. rao 




(3.25) 


The angle 5^ is a design factor to yield the optimum pitch angle 
gain to prevent blade motion instability. Blade pitch, control 
consists of collective pitch control and cyclic pitch control, 

Oc)h a 0.Ct) + ©KiCtJ cW* + ©,.« , (3.26) 

Therefore, the perturbation quantity associated with the angle of 
attack of the blade is simply derived from Eq. 3.7 as 


Sof * 0»f 4* 9^ 


(3.27) 
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3.2 Wing Aerodynamic Forces 


Wing aerodynamic forces are derived by the blade aero- 
dynamic perturbation method. Aerodynamic coefficients are ob- 
tained from strip theory (Ref. 15) . 

The wing element lift, drag and pitching moment (Fig. 7) 
per unit span length can be written as 

a L = i P ui a w 0 f w c w 

^ = if U®- Cd4„ Gw 

U w C Mj)f W ¥ ) C T ( 3 . 28 ) 

These values are transformed to the x,y,z coordinate system: 

F*x - — dO-*#** k. Ar . . 

* -j/’C-w^wr Uir<*w COO^p- - pCwCvow '***«• Ar 

V-*: st _ JD corO <£p. — alL. 

— — ^CwrCft^ Up- C oo<j>^ 

My = -+dL «Cnr 

* "A C ^*0 + ^Jw- 

+ ~^jT f cj- c" (X-w Uw- &W (3.29) 

where 


u; 

Oiff 

-bx^ 4 w 


(V + 1 ^ + U^) a +(u v + Ue,)" 

3nr + pw — 4w 

Uw 4 U< 

V + 4 u& 


( 3 . 30 ) 
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and e is the distance (nondimensionalized by wing chord) between 
the elastic axis and the quarter-chord line (positive if the 
quarter-chord line is ahead of the elastic axis) . 

Hereafter the small angle assumption may be used for the 

wing aerodynamic forces because of the small angle of attack and 

angle of twist. Then it is approximated that sin <}> = 4> , 

w w 

cos <j> w = 1, etc. 

Perturbation velocities u , w . p . u_ and w„ are chosen to 

W W w b b 

express the aerodynamic forces. The perturbation equations are 
then given by 



This study is concerned only with deviation from the equi- 
librium state. Then the steady state aerodynamic forces (F x ) Q , 

(F ) and (M ) may be dropped from the equations without incon- 
z u y u 

sistency although the blade steady state aerodynamic forces (P ) 

z o 

and (P Q ) cannot be neglected since they influence the pylon 
y o 

motion. 

Applying Eq. 3.31 to Eg. 3.29, the following expressions 
for the wing aerodynamic forces due to wing motion are obtained 
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ipc.Va, 

fc,Va 
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M 

0 
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^ />C, V 2 <*w 


Ul 



/’ft* V d-w-Qfw 




(3.32a) 
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— J PC-i* V" OCwb 
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[M 

-/UY Cta v 

< 

W„ 
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P.J 
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^ PCwV ^Tfo 


f H 

o 

4 


— Pc-w *V Cpo^. 


. ** 


(3.32b) 
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M r 

/* V ( Cki* + 

i 

• 

f w ( 
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Pw j 
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u. 

+ 

o 




^■/'cjv i (c«. + a.e) j 


[ivj 


-zPCv v Cc*„ + Offi) 

T 

hi 

o 

- 

M 

/»<£ V( C*. + C* w <*r. + a*«*. e j 

\ 



(3.32c) 

where the angle of attack in equilibrium state is given as 

(y Wo = 0 W (3 ' 33) 
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SECTION 4 


MODAL ANALYSES 


4.1 Variational Functional 

In order to apply Galerkin's method the variational func 
tional is derived from Eqs. 2.4, 2.5, 2.6, 2.7, 2.8, 2.10, 2.11 
2.12, 2.13, and 2.14. 

TT = - L [ 

- & [ 

- £ [ T C( <EI) c V-* 4r ] H 

-u±ifrm+rmn 

+£[ + Uk ,+ flL 

+ £ £-j*J w ^ri» vV- - C°* ) ( K + £ ^>) 

+ >W* 

- (ur„ cWj 4 - 2fLuT H 

+ i^{A‘ + A\^ + 2 a4^| 
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~ x +( e It)b '*+~- 1 Qjj ( By?) 

- Y ^|(HU c -(EI^^-«* <*«G.} 4f 4$**# 

4H^ 

+ u £r < ^^- 

+tC^p^ + ps* 

* 9 

+ "T Mg ^ + -- M £ +i I W$ 

+ A T ?p^P 4 4 1f r v£ 

+ H[f (fi «H. 

+ Z [ j/- Pe *V P*)^ V 

+ z[(V v 'pPa'^4- + ^^“ w ^-' + &)^] 

+ z n*{ - *■ (4«s. - ^ p *h *■ 

■f /04^«.i M + Bt 

-P e uT H «ro<J H | «*/! vL 


27 



♦*[Ci« 4-V , pPg - ^ Es 

+ 5.[iV*H.* 

+ j o Fx ^ d f 

+ L p * ^ ty 

+ C M ^ a » 


(4.1) 


, 4.2 Free Vibration of the Blade and Wing 

In preparation for modal analysis of the basic equations, 
the free vibration of the rotating blade and the cantilever wing 
will be considered. It is very advantageous to use the natural 
frequency of the free vibrating beam to represent the quite com- 
plicated stiffness terms of the beam. 

The variational functional for free vibration of the ro- 
tating blade is derived from Eqs. 4.1, neglecting aerodynamic 
forces, pylon and wing motion. 
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- lJ 2 ( {( e I)c cjw10 » ("a^) 


n r /* 

+ Z, -M mCu^ + vr^ 

K-l k ^ 





(4.2) 


Assume a series solution in terms of normal modes which express the 
coupled motion of inplane and out-of-plane motion 



Substituting Eq. 4.3 into Eq. 4.2 and applying Lagrange's 
equations as follows: 


dt v >jj 


AH* -o 


*% 


Kt 


yields 



<ei) c +(eij, ^*e.} (-§pr)(^$) 


4 {<er) t <^9. + (ei)» ^ e.j (-^s) (-^zr) 


(4.4) 
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+ "5 (CEI) c - CEI).) A^.da (-^ 3 -^ + -^-$) 


+ t -7F-^ +t#-$ +n.*^ViVj]ay- 
» r R 

= Aj j vn(WiWj+ViX|)^ 


(4.5) 


where X . denotes the jth rotating undamped natural frequency of the 
blade. Simplification of Eq. 4.5 will result from use of the ortho- 
gonality condition (Ref. 15) of the normal modes which is expressed 
as 


£ 

f + ViVj)drs 

*4» 


f J*M (V? + VL-)<*r 

* Q 

O 


j*L <4.6) 


The amplitude of the mode shape is normalized to unity as 


& 

[ W C Wl + V*) dr = I 


(4.7) 


Hence, Eq. 4.4 results in 

£[{ (ei) - + (ei) * 

+ a ^ -^9 


30 



j * l 

( 4 . 8 ) 

In the similar way the free vibration of the wing will be analyzed, 
treating the rotor and pylon as lumped masses at the wing tip. The 
wing deflection can be expressed as 

f*j(3)| 

Z I 5j(K) ) (4 - 9) 

a ( W 

Then, pylon motion is described by 

I V« 

'Vr 

* 

l Mr 

= [ "^] ( a ^j l4 * 10> 
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where 


M - 


H= 


V« H • • • H 
4^ * ♦ * 4L. 

<*a <*& Ity 

<&. * . - <f H 

_ dffx d3n 

*3 a 3 * * * ~ d* j 

a,c+j 

Gift) 


{ 4 . 11 ) 


^ cu.fr) 

The wing coupled modes are normalized as 

£ («-<*;-* Jj) + + 

+ K<*J + *J> + if, (■§■)' + V *J + V(t$? 

+ t'(-^' + 3 ;]^ ^ 


( 4 . 12 ) 
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Hence, 



where is the jth wing undamped ..natural frequency. 

4.3 The Variational Functional Described in Modal Form 

The simplified variational functional expressed in terms of 
natural frequencies is 
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- *■ If ($*i JU 

' 1 ( S w j H )(? <#> £ aj ] *"] 

- T 2 ( A X 1 al) * %, a? + i»m 

+i{C«?’S*s>*l 

+ IIC p »<*w< r I 

+ W w , + W v|w (4.14) 


where denotes the work on the wing done by rotor aerodynamic 
forces and W the work on the wing done by wing aerodynamic 

ww 

forces. Details are discussed in Subsection 4. 5* 

The principle of virtual work is applied to Eq. 4.14 to 
obtain the modal equation for the blade and wing* 

For the blade the modal equation can be expressed as 


U + x j h + vj * -v.j(S +Av p )^^ 
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4- -w 0 j Ik. - Voj «*> <l K + Wjj vj 

~ Wjj v) p coo4 K 4 2-ft-Wj 

+ 211 ^ 1 ) 4 , 

= % + p *j 

where 

w y = R'^fr^^Wj^] 

( 1 = 0 , 1 ) 

v a- - R'-^f^v^r] 

u-©, 0 

&J - 

Prj 53 | Pa -vlj 4^" 

For the wing modal equation can be expressed as 


(4.15) 


(4.16) 


(4.17) 
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‘di + A\ a; + 1. - ^)l]^ 

+ 1, s[- v 0 - Vi 


+ w * * v s 


- w u ( ?«j ^-4.) (^-) 3 _- v a - ( % .J W.) (*.J^ 

- 2-Q.Wy "^(g^ -2ilWj OW^ 

= ii 4 <4.18> 

where f i denotes the loading due . to the rotor aerodynamic forces 
derived from the work W_ ro ,. and g. the loading due to wing aero- 

Wr» J, 

dvnamic forces derived from the work W_ — , 

«VK 

4 . 4 Rotor hon-Kotating Coordinate System 


4.4.1 Rotor Non-Rotating Coordinate System 


The equations of motion 4,15 and 4.18 of the blade and wing 
have periodic coefficients because of the appearance of sintj; n and 
cos^. To avoid difficulties due to these periodic coefficients, 
the Fourier-type coordinate system will be chosen to describe the 
motion in the non-rotating system instead of the rotating system 
(Ref. 16). 

(a) For the hingeless rotor, the Fourier coordinate system, 
is described in general as 

%«j ‘ Qj. + S [ ©jc<»> “"MJ +<%,(-> 

+ C-0 H <4.19) 


36 



where Q. is the collective motion, Q. , » and Q . , . are cyclic 

jo ' ^jc(m) js(m) 1 

motion of mth order, defining tilting or warping of the rotor 
plane and Q,* is differential collective motion (only occurring 
for even N-bladed rotors) . 

Since many recent proprotor designs have three-bladed ro- 
tors, this study is concentrated on the odd H-bladed rotor model. 

Therefore, the last term Q, of Eq, 4.19, the differential col- 
• • ^ ^ 

lectrve motion, will be dropped hereafter. In addition to this, 
the higher order terms Q jc(2) , Q js(2) , ...Qj c(m) , Q js(m) are not 
coupled with pylon motion. They represent internal rotor motion 
and usually make only a small contribution to rotor motion. 
Therefore they are also truncated in this study. Hence, 

% k j » Qjo + Qjc c ^4k -h Qj, A4, (4.20) 

Bv inverting the system of Eg. 4.20. the new decrees of 
freedom of the non- rotating system can be obtained: 



The same approach may be applied for the rotor generalised force 
due to periodic aerodynamic forces. 

( Pej + PejOo * ” 5 + 8cj)' 

(Poj + Pyj)s * % X. (£oj •+ P*j) (4.22) 



(b) Tor the gimballed rotor, special care must be taken 
since it has both collective mode shapes (symmetrical modes with 
clamped boundary condition at the root) and cyclic mode shapes 
( antisymmetrical mode shapes with hinged boundary condition), 
while in the hinge less rotor case both collective mode shapes and 
cyclic mode shapes are the same. Therefore, collective motion is 
described with collective mode shapes IJ?, 


,° .rO an( j C y C n c motion 


c c 

with cyclic mode shapes W ^ , V ^ . 



( 4 . 23 ) 


This expression is consistent with the derivation of equations of 
motion up to now, if mode shapes W^, of the coefficients of the 
blade collective motion are replaced by the collective mode shapes 
W °., V? in Eq. 4.15 and 4.18, and those of cyclic motion are re- 
placed by the cyclic mode shapes W?, V^. Therefore, the derivation 
of equations of motion hereafter is based on the hingeless rotor 
case, and those equations can be applied to the gimballed rotor 
case as stated before without confusion. The coefficients asso- 
ciated with collective mode shapes are denoted by superscript 0 
and those with cyclic mode shapes by superscript c. 


The truncation of higher order terms and of the differential 
collective motion in Eq. 4.19 is also applied to the gimballed ro- 
tor case by the same reasoning as for the hingeless rotor case. 


The orthogonality condition for the gimballed rotor case is 
different from the hingeless orthogonality condition. The entire 
rotor system including periodic functions sin^ n and cos<|> n must be 
considered. For example, the orthogonality condition between the 
collective mode and the cyclic sine mode is described as 
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Since 
n=l , . 

4.4.2 


gives 


where 


2. ( * [ wf wj + vj vf 

- { f* « (wfvf 4 

= O (4.24) 


^ is the phase angle defined in Eg. 2.3, the summation over 

. . N of sinij; becomes zero, 
n 

Pitch-Flap Coupling in the Non-Rotating Coordinate System 
Pitch-Flap Coupling is expressed in Eq. 3.25 as 

v---!£| (3 - 25 > 

Wad 

The substitution of mode shapes into this expression 


0 PF s, 2L[- 

= Kp* [ Qco Qi^c***^* + Qi* (4. 


25) 



(4.26) 
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For the hingeless rotor case considered here, K pi is always 
zero because the slopes of blade mode shapes at the root are zero, 
due to the clamped boundary condition. 

For the gimballed rotor case, Eq. 4.25 is rewritten as 


where 



( 4 . 27 ) 


K 




( 4 . 28 ) 


From the definition of the collective mode shapes, K ^ is 
zero. Therefore, for the gimballed rotor case the pitch-flap 
coupling contributes only to the cyclic motion of the rotor, not 
to the collective motion. 

4.5 Aerodynamic Forces Described in Modal Form 
4.5.1 Blade Aerodynamic Forces 

From Eq. 4.17 blade aerodynamic force is described as 

r* 

J?*j = \ (PeVj ( 4 . 29 ) 

The application of Eq. 4.22, accompanied by the substitu- 
tion Eq. 4.21 into Eq. 3.13 yields 
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(Pej+Pcj). = T J I‘ n t G J L 0; " + ^ 

+ ^Jtj ~ + A <£ ay + Svj -5 k*. ^-o 

+ $V/ &] (4 * 30) 

(P«j *f P*j) t = J- y x B rf[ Q; c 4 ja X, <$ji Q: s 

~ ^ + * ^ x j °V ~ G*rj V p 

4 A $xj --■ -? 4 ^ Q.'c 

+ «**] M.31) 


(% 


^ S’* ^‘ e ~ “j^r ” & 5x f * v£ 


O 

5 ^p 


+ A <Sr,' Wf> + ^ ijf 

+ 'Syj- 9,jj 


(4.32) 
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where 


1 

.[' ( p 2SL2S 
4 (P ®' Q* 

+ 1 

- >£i2£ + 

’*1 R i 


. r- X;^c 

+ F ia ~^r- 


* 


» 

/~N 

fl 

p|<5? 

-4- 




r f vj 

~J + 

■ / o 


•, -^} scrfac 

<*3B.j 

‘l'^R + 

FV a -^)c4^ 


S *J -( (P»j-r + Fi^)xdz 
G Vj - Jf ( Pbj-j^- + -7^) dx 


*»• “ &•»/ [ i > Ib 0%) J 

*' = ($!/[* 

*•» ' (^./[VlB-^J 

Rj • (#3./ [w^n 

Fm " ■ ft "/ [i:«b(4-)1 

^ “ (^./[iu^] 

FlJ * (4fH/[in. -^j 

F « - 


(4.33) 


(4.34) 
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4.5.2 Wing-Loading Due to Blade Aerodynamic Forces 

The virtual work done by wing loading due to blade aero- 
dynamic forces is described as 

* | %r x %r t x 


” fo ^ Be 

Wp 4-v^ ^ ^ 

--&.(& + V ^> Y '-P e + ^ ^ <=«4 n - ^ i*r„ c*s c l^ 


|f& 


C4.35) 


The application of the blade mode shapes and non-rotating coordi- 
nate system gives 0 


Sw * 

WB 


* u 


( Sr v6 

T 

* 

r Vr 


/Vr' 

S Vr 
S v s 1 

k 


Vr 

v 3 

' 

^ +k[a3]i 

v« l 

Vu ( 

Sv r 

l S*V j 


i ! 

t l/ r J 


V 

vV 

k Vy J 


+ ?«[AaJ|o ; ] +s*[a?J{o 1 | 

IVv 

1 w 

+ *[*£*]{ *V > 

1 60 I 

I I 

V J 


{4.36) 
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where k is an integral operator defined as 


K 



dm 


and integrand matrices are 



_ 5 L J*i 1 ° j 

-ixr„ J _o_ T 2 «•> 

-^•i.Fei j^o Jji^x(-T : ®x+ F’j,') j^ 
o j xF et | o | 



| 4 -lFej_ _ j _ G_ 

1 O 

| 4:^ x (^a -^ii) I O 

I c> 

I _ o 


o| o 

o_ j__ 

o I o 

— I- - 

O I 0 


, 0 | -iFe,+Ti 0 jo 

_J_ 6 I _ _o_ _ i o 

1 * l_ Jxf*, «5 *J[o 

J -T*** _ j_ 'Uyz a ^)_l O 

10 I o \ 0 . 


(4.37) 


(4.38a) 


(4.38b) 
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[4j= 


Vi r , VCt- 

o 

o 


+ -§ fF «0 


4 


o 

o 




“ C % F *l ^R ^a) I “1 (§ Pft, +^‘ ^ 

r -_- 


(4.38c) 


M5- 


° ! 

| 

i \ ** Pe 3 

*¥1^ j O 

j o 

I ^[(#w-$h 

1 4&(^f^ + i£ ; Rh) 

o ; -rLx(^r*, + ^>F^)j 


] “i -tKpl Fej 

1 

1 ^ ^ FoJ 

0 | 

1 

1 

1 " S-X ( l^T^) 

1 

l_ 

] -i ^^3 


1 ° 


K*] 


(4.38d) 


A r- 

2 nil 


O | o o 

. l-?5^ JjSL 


I- 


t o I 

I I 

I 0 1 o 

f- 


1 45^1 1 __°_ _ j _o _ i -^tR* i ^*xF«3 

4 **Ll 4 l*J* I <> J 

o I ° I J xFtki ] o o 


(4.38e) 
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4.5.3 Wing Aerodynamic Forces 

The virtual work done by wing aerodynamic forces is ex 
pressed by 


= | S S wjjr 





( 4 . 39 ) 


The displacements u^, w w * p w are described in terms of the mode 
shapes as 


< 




where 



awn 


>. r..\ 

• • • 

5, on 


■ • • ^hQA 

$ca) 


• • * 


( 4 . 40 ) 


( 4 . 41 ) 


Then, the substitution of Eg. 4.40 into 4.39 gives 



( 4 . 42 ) 
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where 


[ <d £ [sf [a ,1 [s]^ 

[Af J'TKrC [5lT ^ d 3 


The matrices A , A and A_ are derived from Eq. 3.32 as follows: 

W W (j 


[A,] 


-±pC m Va 9 (l+ j pCmVGvt vw o 'O 

I j — 

j _jo G,V I ° 

-^pci VCC^-H^a ) ! pd VCC,. + C* I 

1 I ° 

i +a-w<vw.jc; > 


(4.44a) 




0 I -LpCmV'a* 

I 



0 { Cw V (G* ^ 

I 


(4.44b) 
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j o ; PtV*Q,< 


- 4^-f 1 d w C/*„ j q I - P Cw"V^Ca »»■ I 

-^/>cJV(Gn,+ fl*?) I ° | PGi v\a..+c.^, j 

| I + cw*war) ! 

i * i 

[ o ! o ! o 

i L i 

■ I I 

[ o ; o ] o 

I - ' 1 1 

! o ! o ! o 

, i : 




Final Equations of Motion 


In this part the final equations of motion will be derived. 
They will be expressed in general format, so that any number of 
mode shades both for the blade and wing can be chosen to describe 
the rotor or wing motion as accurately as required. 


4.6.1 Equations of Motion in General Format 


The substitution of Eqs. 4.20 and Eq. 4.10 into Eq. 4.15 
gives the jth blade modal equations of motion. Note that time 
derivatives are normalized by rotational speed as in Eq. 3.24. 


For the jth blade mode 


oo © . . 

{Qj} 4 h f 4 [ M ®]j \ Q J| 

4[Mj][T]{aV 5 ’[ c j]fT]^i < 4 . 45 , 
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where 





(4.46a) 


(4.46b) 


(4.46c) 


(4.46d) 


(4.46e) 


(4.47) 
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where 



50 



Loading matrix^ f} for the wing consists of aerodynamic 
forces f^ due to the rotor, and aerodynamic forces g^ due to wing 
motion itself. 



Similarly, aerodynamic forces are derived. 

For rotor aerodynamic forces 

(&j+P*j ) 0 

= 1? I 4 / 

l CPof f Pzj) 3 j 

■?f41j (i| 1“.-) 

*W!j[Tl|a|.[A4[T](a} 

* |*| 

where 
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(4.52b) 



(4.52c) 


(4 . 52d) 
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0 . A Gfj o o 

0 ° % 0 
o o o 0 

(4.52e> 


The excitation vector (e) is expressed as 



m 


V 






(4.53) 


For wing aerodynamic forces 
From Eg. 4.50 

J f| = { £ + g| (4.54) 

Loading to the wing due to rotor aerodynamic forces f i cor- 
responding to is expressed from Eq. 4.36 as 
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{-f}- <[Tf[A“][T]{a} + K [t] T [ Aa'] [t] ( a J 

C. 

+ ? «M T [/&]{&;} 

-<W T [ArjH 

Loading due to wing aerodynamic forces g^ is described from Eq. 
4.42 as 

{$j =[a“] (a} + + [A?]H ( 4 . 56 ) 


4.6.2 Equations of Motion for Nine Degrees of Freedom 

For a typical case, nine degrees of freedom are considered 
here. These degrees are rotor flapping motion (collective flap- 
ping motion Q 10 » cyclic flapping motion Q lc and Q lg ) , rotor lag- 
ging motion (collective lagging motion §20' cyclic lagging motion 
Q and Q 0 ) , wing vertical bending motion , wing chordwise 
bending motion a 2 and wing torsion a^. 

These mode shapes are natural modes, therefore they are 
coupled . 

The final equations in matrix form are 

[a] {*{ + [B] <*J + [c] |*j - [D] {«J ««•«> 


54 



where 



( 




O.0 \ 

i 

f Ufi. \ 

a* 


V 

Jik 

Q»s 

I 

V 

Oio 

\ w= 

ur 6 

v 

Qic 


9. ' 



9ic 

a, 


e,J 

ax J 



a 3 / 




The coefficient matrices [A], [B] and [C] are 9 
and [D] is a 9 x 6 matrix, as given feelow. 
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(4.58) 


9 matrices. 


(4.59a) 
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with Uncoupled Modes 


In order to have a better understanding of tilting prop- 
rotor dynamics, a simple special case was considered. This case 
is based on the assumption that the mode shapes are uncoupled. 
Therefore, rotor motion is expressed as 


(**1- | *M(q..+ q, 

[\r H ) 0 ) 


4 « + 


+ { v j (Qao + ^4.*) (4.60) 
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New symbols will be defined below, but all these symbols 
can be derived by applying the uncoupled mode shapes to the general 
equations. 


Pylon Motions 



> i 

at 

ii 

:P 

'Wt » 

s* 

at 

y- 1 - 

< = 

dXx 

at 

L 

I = 

<*>> 

at 

3-L 

w; = 


at 

* = «- 


(4.63) 


Therefore the matrix [t] ih Eq. 4.10 is described as 


fTj = 


Via 

R 

O 

0 

O 

Hi 

Q 

o 

O 


o 

0 

0 

$ 

-W* 

o 

o 


Aerodynamic Forces 

(a) Rotor Aerodynamic Forces 

Pjf — i*!. ^ ^r.) A*. 

Lis =-T > $: CPbJ 

Ljj --id, ^ ( F.,) V_ £ X 


(4.64) 


(4.65a) 
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(4.65b) 


(4.65c) 


(4.65d) 


'-pp = ~ ^ ^ fc CR*i)^ L ctt. 
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where in general F,L expresses the forces corresponding to the 
direction of flapping motion and lagging motion respectively. Sub 
scripts 8 and £ describe the mode shape contributions to the aero- 
dynamic forces? T and P show the derivatives for the perturbation 
expression with respect to inplane velocity and out-of-plane 
velocity U ; superscript + expresses the moment at the blade root. 


cr < ,= 4rr»^CF ao )«ix 

%^nr e [; 
u--in 6 

L P - x 
L* s ■j f I® 


Fp + « 4; V I® ^ ( Ft») 3^1 

i-T — yi* 

Le = i*** CFes)*^ (4.66) 


c a„ = -M c: 

Ca„ = -M [ -Up-«i -Si _ T5* 

c a.i “ 
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c«„ = - a t + i ( S'lt -I- J(m3] 

Ce M ■ -N [ih^-F?)w;i] 

S- - u [4? 1 -r5 : ^J 
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Co,» - -N [ ic F. - 1 l t ) (irif ] 
Ca„- -nUt*? +Lt)^€j 
C« M - - N t"S F J 

C a44 . -Mt 


(b) Wing Aerodynamic Forces 
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(4.67b) 


(4.68a) 
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c ik C a* »4j ix*% 


^" a <» 25 

C *3i = ” £* 5„ 


C *" s "c? L A «« 5f ' < *^ 

^ - i*. 1 ****“*3 

** jjf ^ A§n 5-*-^tJ 

^ " fi5 C^**^*^ (4.68c) 

Uii - -5? S. ***+**!} 
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where A. , A , and A„ are components of matrix [A.], [A 1 , 
13 iJ 13 

and [Ag] respectively in Eq. 4.44. 

Pitch-Flap Coupling 

The pitch-flap coupling coefficient is expressed as 



<AvJi 





(4.69) 


The rotor collective mode shape has zero slope at the root. There 
fore, pitch-flap coupling appears only in cyclic motion. 

Natural frequencies 

X : blade flapping natural frequency 

P 

X^ j blade lagging natural frequency 

: wing vertical bending natural frequency 

: wing chordwise bending natural frequency 

A 3 : wing torsion natural frequency 

4.7 Equations for Gimballed Rotor 

As discussed in Subsection 4.4, equations of motion (Eqs. 
4.45 and 4.48) are applicable to the gimballed rotor. The signi- 
ficant difference between the gimballed rotor and the hingeless 
rotor lies in the application of the blade mode shapes to the 
rotor motion. The gimballed rotor has independent collective 
mode shapes for the collective motion and cyclic mode shapes for 
the cyclic motion. Therefore, the above distinction is found in 
the following matrices: 

For blade equations 
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(4.70) 



(4.71a) 


(4.71b) 
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(4 , 71d) 
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SECTION 5 


NATURAL FREQUENCIES AND MODE SHAPES 
OF THE ROTOR AND WING 


5.1 Rigid Rotor 

In order to obtain the natural frequencies and correspond- 
ing mode shapes of the free vibration of the rotating blade, the 
equations of motion are derived neglecting aerodynamic forces and 
pylon motions in Eqs. 2.3 and 2.4. 

I* [(Me + (El). c**'Q*} 

+ [J CE-£) c - Ce t) B j 0 S aw 6a 

-§r( T #) +*v^=° 

|L £ j CEI) c <***8 e + tel) a ©a) “1^1 

+ .|L [ j CEr) c - Cei)a} ] 

- ■fr( T 4^') 4 ~ ^ ° (5.1) 

These equations of motion are rewritten as a usual eigenvalue prob- 
lem to find eigenvalues and/or eigenvectors. 

Boundary conditions for the rigid rotor are 
at the root: 

= O 
tr; 



(5.2) 
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at the tip; 


i** =o 

=o 

ar 4 


ar* 

._ 


arA 


= o 


(5.3) 


As a typical hingeless rotor, the Boeing rotor is chosen in this 
study. Mass, stiffness distributions and built-in angle of twist 
are shown in Pig. 8. 


Nonrotating natural frequencies of the blade are shown in 
Fig. 9. Almost no influence of the collective pitch change on 
the natural frequency is found in the case of the rigid rotor. 


In Figs. 10a and lob, the first and second natural fre- 
quencies of the rotating blade are shown versus rotor rotational 
speed and inflow ratio. Inflow ratio is related to the collective 
pitch as described in Eq. 3.9. At normal rotational speed (in 
this case ft = 386 R.P.M.) the natural frequency variation due to 
collective pitch change is shown in Fig. 11. Typical mode shapes 
up to the 4th mode are shown in Figs. 12a through 12d. Blade mode 
shapes are normalized with respect to the blade radius at the 
point of maximum deflection in either out-of-plane or inplane 
bending. 

5.2 Gimballed Rotor 

In the case of the gimballed rotor, rotor motion is ex- 
pressed by collective modes and cyclic modes (in Ref. 17 and 18). 
Hence, based on Eq. 5.1, the eigenvalue problem should be solved 
with boundary conditions for collective modes and cyclic modes 
for the powered flight, respectively. 

Boundary conditions for collective modes are defined to 
yield symmetrical modes for flapping and lagging. Therefore, 


at the root : 
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(5.5) 


__ 


= o 


a*-* “ ar* 

Boundary conditions for cyclic modes which consist of anti- 
symmetrical modes for flapping" motion and symmetrical modes for 
lagging motion are expressed as: 


at the root: 


f, 
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=o 
av*- u 


^ - o 

=o 

(5.6) 


air 


at the tip: 



“3?r =» 0 

Ju L 


0 

It 

$ 

- ° 

(5.7) 


^2 ^ typical case , the Bell design is considered here# Mass; 
stiffness distributions and built-in angle of twist are shown m 
Fig ♦ 8 ♦ 

Blade nonrotating natural frequencies for the gimballed 
rotor are shown in Fig. 13. In Fig. 13b, the first natural fre- 
quency for rigid body flapping motion has non— zero natural fre- 
quency due to the rubber hub spring, which is intended to increase 
control power and damping. 


In Figs. 14 and 15 it is shown that collective pitch change 
has a large influence on the natural frequency variation except for 
the first cyclic mode natural frequency (rigid body mode) . Mode 
shapes for the gimballed rotor are shown in Figs. 16 and 17. The 
normalization system is the same as that of the rigid rotor. 



5.3 Wing 

The equations of motion for free vibration of the wing are 

the same as Eq. 2.5, 2.6, and 2.7 when the aerodynamic forces F x , 

F . and M are eliminated from these equations. Boundary condi- 
z' y 

tions are as follows: 


(a) at the root : 

— 0 
uJiv = 0 

P«r = O 

(b) at the wing tip: 
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Xpp^p + ] _ u +^Ma^.r,c 

+ MMe^p +*^18^=0 (5.9d) 

V £ ■ [| (F1 ^c + CEr,)a cao ^}4§^ fsU 

- [((EU-CEiJ. -|&] ^ 

+ NXiH. -0 (5.9e) 


The relationship between pylon motions and wing deflections is 
defined in Eg, 2*8. The rotor and the pylon are treated as 
lumped masses at the wing tip in these equations; as in the actual 
structural dynamic test, the proprotor blades will be removed and 
replaced by equivalent weights. Note that in Eq. 5*9 the blade is 
treated as a lumped mass which has equivalent weight and equivalent 
blade flapping inertia Ig. This flapping inertia leads to lower 
wing frequencies, especially in torsion, than when the blade is 
treated as a lumped mass with equivalent mass and without equi- 
valent blade flapping inertia. 

The same wing is used for the Boeing and the Bell design; 
mass and stiffness properties are shown in Fig. 18 ♦ The differ- 
ences appear in the pylon and blade mass properties shown in Table 

1. No built-in angle of twist, dihedral angle, and sweep angle are 

considered here. 

Natural frequencies of both cases are listed in Table 1, and 

mode shapes are described in Fig. 19 for the Boeing case and in 
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Fig. 20 for the Bell case. Wing mode shapes are normalized with 
respect to the wing semispan at the point of maximum deflection 
in either vertical bending or in chordwise bending. However, 
the mode shape corresponding to the third natural frequency of 
the wing, which is predominantly one of wing torsional deflec- 
tion, is normalized at the point of maximum torsional deflection. 
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SECTION 6 


PROPROTOR DYNAMIC CHARACTERISTICS 
6 . 1 Introduction 

In this section the eigenvalues and frequency response of 
the system will be discussed. The results are compared with those 
of Ref. 13. The analyses were conducted for two types of prop- 
rotor design: the hingeless rotor Boeing design, and the gimballed 
rotor Bell design. The case selected to be investigated is cruis- 
ing flight in the airplane mode for each design. The data for the 
calculation are listed in Table 1. One should notice that the 
equations of motion shown in Eg. 4.57 were derived based on mass- 
normalized coupled modes. However, for the sake of aiding physi- 
cal understanding of the frequency response and eigenvalue analy- 
sis, those results are based on physical mode shapes, as presented 
in Figs. 12, 16, 17, 19, and 20. 

The primary differences between the Ref. 13 analysis and 
this report are tabulated below. 

This Report Ref. 13 

Natural Frequencies Calculated Experimental 

Data Data 

Coupled Modes Uncoupled Modes 

Coupled Modes Uncoupled Elastic 

Modes for Inertia 
Terms ; Rigid-Body 
Mode for Aero- 
dynamic Terms 

Structural Damping No Yes 


Mode Shapes 
Wing 
Rotor 
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The case studies done in the present study are tabulated 

below: 


Degree-of- 

Freedom 

Mode Shapes 
Blade 

Wing 

Analysis 


9 Degree-of-Freedom 18 Degree-of-Freedom 


Rigid Body 
Coupled 


Uncoupled 

Elastic 

Coupled 


Coupled 

Coupled 


Coupled 


Coupled 


Eigenvalue Eigenvalue Eigenvalue Eigenvalue 

and Eigen- 
vector 


Frequency 
Response to 

V V G' W G' 


The general symbols Qj C r Qj C » ® js ' a i ^ or b ^ ade and wing 

motions are convenient for the theoretical derivation of the equa- 
tions , but from the physical understanding aspect, it would be 
better to choose other symbols. For this section, except for the 
eigenvalue analysis, the following system will be used to express 
blade and wing motion for the 9 degree-of- freedom system: flapping 
and lagging motions which are designated as 6 and ? , respectively, 
and wing motions are , q 2 , and p, which express the vertical 
bending, chordwise bending, and torsion, respectively. 

The correspondence between these sets of symbols is tabu- 
lated below: 
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New Symbols 


General Symbols 
Bell Design Boeing Design 


e n 

Predominant 

°2a 

°2o 

O 

Flapping 

®lc 

Motion 

Q lc 

Q 2c 

6 ls 


Q ls 

Q 2s 


Predominant 

Q 10 

Q 10 

o 

Lagging 

5 lc 

Motion 

q 2c 

Q lc 

C ls 


q 2s 

Q ls 

*1 

Wing Vertical 
Bending 

a i 

a l 

<*2 

Wing Chordwise 
Bending 

a 2 

a 2 

P 

Wing Torsion 

a 3 

a 3 


6.2 Eigenvalues and Eigenvectors 

The system stability characteristics are shown in Figs. 21 
through 24. 

6.2.1 Addition of Higher Mode Degrees of Freedom 

The eigenvalue variation with number of degrees-of- freedom 
considered is shown in Fig. 21. The 9 degrees-of- freedom include 

Q 10' Q lc' Q ls' Q 20' Q 2c' Q 2s' a l' a 2 ’ a 3 (B o' 3 lc' 6 ls' C o' 

^ q^, and P) corresponding to the blade’s first two 

natural frequencies and the wing’s first three natural frequencies. 
This system has 9 eigenvalues and 9 corresponding eigenvectors. 

The 9 eigenvalues are denoted by their frequencies and the partici- 
pation of the degrees-of- freedom in the eigenvectors is indicated 
below: 
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Present Ref. 13 


Eigenvalue Symbol Participation 

Bell 

Boeing 


Low frequency 

6-1 

C-l 


Collective 

5 

? 


High frequency 

6+1 

C+l 


Low frequency Q 2 

c-i 

6-1 


Collective Q 2 

6 

6 

Q 2 

High frequency 

C+l 

6+1 

a i 

The first wing 


q l 


frequency 



a- 

The second wing 

q 2 

q 2 


frequency 



a~ 

The third wing 

P 

P 

j 

frequency 



In the 

18 degree-of -freedom system two 

higher 

elastic 

modes for the 

blade and three higher elastic 

modes for the wing 


are added to the 9 degree-of- freedom system. The same designa- 
tion system as used for the 9 degree-of-f reedom is employed for 
the higher eigenvalues. The results of the 18 degree-of-f reedom 
analysis show that the addition of more degrees of freedom does 
not substantially influence the basic system eigenvalues since 
the added degrees of freedom have large natural frequencies in 
comparison with the original values (see Table 1) . The 18 de- 
gree-of- freedom system eigenvalue locations are shown in Fig. 22. 
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6.2.2 The Sensitivity to Mode Shape 


Eigenvalues and damping ratios (the fraction of the criti- 
cal damping) for three different types of mode shapes for the 
blade are shown in Fig. 23. One is the rigid— body mode (with 
spring restraint at the root) , another is the elastic uncoupled 
mode, and the third is the elastic coupled mode, while the wing 
mode shapes considered are restricted to the elastic coupled mode 
shapes only. 

The results tell that mode-shape types make negligible dif- 
ference in the frequency of the eigenvalues; however, they in- 
fluence the damping a lot. 

In the Bell design, there is almost no difference between 
the rigid mode and the elastic mode in the damping ratio. This is 
due to the similarity between the rigid— body mode and the elastic 
uncoupled mode. In the elastic coupled mode, a slightly higher 
damping is obtained. 

In the Boeing design, each mode-shape type gives a differ- 
ent damping. For the first natural frequency mode of the blade, 
the rigid-body mode calculation is conservative rather than that of 
the elastic coupled mode-shape type, and for the second natural 
frequency mode, it is nonconservative. This is because the first 
mode has both positive out-of-plane and inplane deflection, while 
the second has a positive out-of-plane and a negative inplane 
deflection as shown in Figs. 12 (a) and (b) . The corresponding 
positive out-of-plane deflection of the first mode increases the 
damping, while the associated negative inplane deflection of the 
second mode reduces the damping. This explanation also holds for 
the Bell case. The damping of the second collective elastic coupled 
mode of the Bell rotor is lower than that of the rigid-body mode. 
The second collective mode has a positive out-of-plane deflection 
and a small positive inplane deflection. However, the rotation 
direction of the rotor is different for the Bell and the Boeing 
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designs. Therefore, the positive inplane deflection mode of the 
Bell rotor is physically the same motion as that of the first 
mode of the Boeing rotor. The coordinate system related to this 
is shown in Fig. 5, and the contribution to the damping of the 
coupled mode shape can be seen in the term in Eq. 4.33. 

In conclusion, the use of the elastic coupling mode affects 
the damping significantly. The coupled mode with a combination of 
forward out-of-plane deflection (upward in helicopter mode) and 
inplane deflection opposing the rotor direction of rotation, or 
vice versa, increases the damping over that of the rigid-body- 
mode calculation. The coupled mode which has forward out-of-plane 
deflection and inplane deflection proceeding with rotor direction 
of rotation, or vice versa, decreases the damping. This holds 
only for the first beam bending mode; in other words, the coupled 
mode without nodal points between the root and the tip. 

6.3 Frequency Response 

The frequency response analysis is dealt with in this sub- 
section. These calculations are all based on mass-normalized 
coupled modes for both the blade and the wing. However, for con- 
venience in physical understanding, the results of the frequency 
response calculation are presented in terms of length-normalized 
mode shapes in Figs. 25 through 32. 

6.3.1 Frequency Response to the Gust 

Frequency responses to the vertical gust u„, the lateral 
gust v Q , and the longitudinal gust w G for both the Bell and the 
Boeing designs are shown in Figs. 25 through 30. 

As a whole, the behavior of the frequency response is quite 
similar to that of Ref. 13 in spite of the differences in the theo 
retical model stated in Subsection 6.1. 

The detail characteristics of the frequency response are 
discussed next. 
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6. 3. 1.1 Blade Collective Motion Response { , c o > 

Collective motion responses have a close relationship with 
» the wing motions (vertical bending q^, chordwise bending, q 2 , and 

torsion, p) . They have strong resonances with these modes. In 
the case of and v Q inputs, the static response (u equal to 
zero) and the lower- frequency- range response are negligible. How- 
ever, the response to the longitudinal gust w_ has a significant 
static and lower frequency range response. Resonances of the col- 
lective modes occur in the response to the longitudinal gust input. 
Comparing the Boeing design with the Bell design, the Boeing has 
a larger response in collective responses to each gust input. 

6. 3. 1.2 Blade Cyclic Flapping Motion Response (8^ c and 6^ s > 

For the cyclic inputs (u G and v Q ) there are not strong 
resonances between cyclic flapping motions and the wing vertical 
bending motion (q^) . In the upper frequency range, the high fre- 
quency flapping mode (S+l mode) has a significant resonance ap- 
pearing in the flapping motion. In the lower frequency range, 

„ there is a resonance of the low frequency flapping mode (8-1 mode) 

for the Boeing design and the low frequency lagging mode (c-1 mode) 
for the Bell design. 

For the collective input (w G ) the response of the cyclic 
flapping motion has resonances with the wing motion modes (q^, q 2 , 
and p) in both designs. 

6. 3.1. 3 Blade Cyclic Lagging Motion Response and C^ g ) 

In the Bell design there is an evident resonance in the 
cyclic lagging motion response and S ls ) to the vertical gust 

input (u G ) in proximity to the wing vertical bending mode (q^) . 

To the lateral gust input (v G ) , a low frequency lagging mode 
U-l mode) resonance appears in the lagging motion response. 

In the Boeing design there cannot be seen such obvious 
resonances in the response to the vertical and lateral gust inputs. 


83 


The static response is much larger than that of the Bell design 
due to the soft inplane design (the blade lagging natural frequency 
is less than unity) . 

For the collective input (w^) , the responses of cyclic lag- 
ging motion have resonances with the wing motion modes (q^, q^ , 
and p) in both designs. 

6. 3. 1.4 Wing Motion Response (q^, q 2 , and p) 

The wing motion includes resonances in each wing mode. 
Although the response magnitude is quite large, it is expected to 
become rather small if structural damping is included. 

6.3.2 Frequency Response to Control Pitch Angle 

Frequency responses to longitudinal cyclic pitch input are 
shown in Figs. 31 and 32. 
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SECTION 7 


CONCLUSIONS AND COMMENTS 


7.1 Conclusions 

This study has been devoted to the development and evalu- 
ation of a theoretical model of the proprotor on a cantilever 
wing, operating in normal cruising flight. This theory expresses 
the wing and blade motions in coupled form, and can include any 
number of mode shapes required to express the motions accurately. 
It has been applied to an investigation of the dynamic character- 
istics of the Bell and the Boeing designs. 

Based on the theoretical results included in this study, 
the following conclusions may be stated: 

(a) The choice of mode shape (rigid-body mode or elastic- 
coupled mode) affects the damping significantly. The 
dependency of the damping on the mode shape can be esti- 
mated for the first beam bending mode. The blade inplane 
deflection opposing the rotor direction of rotation, 
accompanied by the forward out-of-plane deflection, in- 
creases the damping, comparing it with the rigid-body 
calculation. The inplane deflection proceeding in the 
rotor direction of rotation decreases the damping. The 
mode shape has little influence on the frequencies of 
the system. 

(b) The results of the frequency response are quite similar 
to those of Ref. 13, in spite of the difference in the 
mode shapes. The amplitude of the response is slightly 
different, since structural damping was not included in 
the present calculation, and the mode shapes used were 
different. 
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(c) The analysis of the eighteen degree-of-freedom system showed 
that the higher- frequency degrees of freedom have small 
influence on the basic degrees of freedom. 

7.2 Suggestions for Future Research 

A direct and useful extension of the present study would be 
the stability analysis of the proprotor aircraft, with respect to 
air resonance and flutter. This theory will be very powerful be- 
cause the eigenvector components can be compared directly without 
any adjustment between components. Mach number effects of the 
blade should be included for the flutter analysis. 

Stresses or bending moments of the wing or blade can be pre- 
dicted from the motions of the wing and blade obtained from this 
analysis. In addition, this analysis may be applied to the de- 
velopment of an automatic control device to alleviate the gust re- 
sponse of the vehicle. 
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TABLE 1 


DESCRIPTION OF THE BELL AND THE BOEING PROPROTOR 
DESIGNS CONSIDERED IN THIS REPORT 


BELL BOEING 

ROTOR 


Type 

gimballed, stiff 
inplane 

cantilever, soft 
inplane 

Number of 
blades, N 

3 


3 

Radius , R 

156 in. 


150 in. 

Chord, Cg 

18.9 in. 


14 in. 

Lock number ,y 

3.83 


4.04 

Solidity, a 

0.089 


0.115 

Pitch/flap 
coupling, 6^ 

-15 deg. 


0 

Collective pitch, 
0 D 

1.25 deg. 


1.0 deg . 

Lift-curve slope, a 

5.7 


5.7 

Drag Coefficient, 
C Do 

0.0065 


0.0065 

Rotor rotation 
direction, 57 

+1 


-1 

Inflow ratio, \ 

0.7 


0.7 

Rotational speed , | ft ] 

458 RPM 
48.9 rad/sec 

386 RPM 
40.4 rad/sec 

Blade Natural Frequencies 



first, X^/|ft| 

1.02/rev 

(7.78Hz) 

0.827/rev. (5.32Hz) 

second, X 2 /I ft] 

1.34/rev 

(10.2Hz) 

1.32/rev (8.49Hz) 

third, X^/l & | 

4.35/rev 

(33.2Hz) 

3.40/rev (21.9Hz) 
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TABLE 

1 CONTINUED 


ROTOR (cont'd) 

BELL 

BOEING 

fourth, A 4 /| & | 

Collective Natural 
Frequency 

10, 1/rev (77.1Hz) 

6.77/rev( 43 

first f X 1 (o, /|n| 

1.31/rev (10.0Hz) 


second, xj°V|ft| 

2.12/rev (16.2Hz) 


third, 

4.93/rev (37.7Hz) 


fourth 

10.6/rev (80.9Hz) 


Blade flapping 
inertia, X fi 

105 slug-ft 2 

150 slug-ft' 

One blade weight. Mg 
WING 

133 lb 

124 lb 

Semispan , L 

200 in. 

200 in. 

Chord, c 
w 

62.2 in. 

62.2 in. 

Mast height, h 

51.3 in. 

51.3 in. 

Sweep 

0 

0 

Dihedral 

0 

0 

Lift-curve 
slope , a w 

5.7 

5.7 

Drag coefficient, 
c 

Dow 

0.004 

0.004 

Moment coefficient C 

mo 

-0.005 

-0.005 

Aerodynamic center, 

5 - V C W 

0.01 

0.01 

Angle of attack, a 

9 wo 

2.0 deg 

2.0 deg 
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TABLE 1 CONCLUDED 



BELL 

BOEING 

WING (cont'd) 

Natural Frequencies 

first, A^y] 

0.347/rev (2. 65Hz) 

0.36 5/rev (2.35Hz) 

second,A2/|ft| 

0.622/rev (4. 75HZ) 

G .653/rev (4 ,20Hz) 

third, A 3 /|fl| 

1.09/rev (8.32Hz) 

1 . 11/re v ( 7 . 14 Hz ) 

fourth,A ( j/|tt| 

2 . 37/rev(18. 1Hz) 

2. 4 7/rev (15. 9Hz) 

fifth, A 5 /|£i| 

3 .76/rev (28. 7Hz) 

3.95/rev(25.4Hz) 

sixth, Ag/jft| 

10 . 6/rev (80 . 9Hz) 

12.5/rev(80.4Hz) 

PYLON 

Weight, M 

P 

1420 lb 

2000 lb 

Yaw inertia, I .... 

164.8 slug-ft 2 

250.0 slug-ft 2 

py 

Pitch inertia, I__ 

190.0 slug-ft 2 

250.0 slug-ft 2 

PP 

Roll inertia, I 

42.4 slug-ft 2 

30.0 slug-ft 2 

FLIGHT CONDITION FOR 

CALCULATIONS, A = 0.7 


Cruising speed, V 

250 kt 

218 kt 

Cruising altitude 

sea level 

sea level 


91 




FIG. 1 CONFIGURATION OF ANALYTICAL MODEL , PROPROTOR ON A CANTILEVER 
WING OPERATING IN THE AIRPLANE MODE 







FIG. 2 PROP ROTOR GEOMETRY AND MOTION WITH GUST VELOCITY DEFINITION 






AFTER 



FIG. 4 WING GEOMETRY 
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X 



FIG. 7 VELOCITIES AT THE WIilG ELEMENT AND RESULTING 
AERODYNAMIC FORCES WHEN LOOKING OUTBOARD 
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(b) Section Flapwise Bending Stiffness Distribution 


FIG. 8 CONTINUED 
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800 


BELL 

BOEING 



FIG. 8 CO 


iNTINUED 
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DEG 


40 U 


BELL 


BOEING 



1.0 


(d) Anglo of Twist 
FIG. 8 CONCLUDED 



FIG. 9 BLADE NATURAL FREQUENCIES, NONROTATING, 
FOR BOEING ROTOR 


E 


3 



NATURAL FREQUENCY PER REVOLUTION 



INFLOW RATIO, V/|fi|R 
(a) 1st Natural Frequency 

FIG. 10 BLADE ROTATING NATURAL FREQUENCIES FOR BOEING 
ROTOR 
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NATURAL FREQUENCY PER REVOLUTION, X 2 /|ft| 


2.0 


1.0 



a = 300 RPM 



= 386 RPM ( NORMAL fi) 



0 I 1 

0 1.0 

INFLOW RATIO , V/ ] H | R 
(b) 2nd Natural Frequency 

FIG. 10 CONCLUDED 


—I — 

2.0 


1U5 



NATURAL FREQUENCY PER REVOLUTION 



FIG. 11 BLADE NATURAL FREQUENCIES FOR BOEING ROTOR AT 
CONSTANT ROTATIONAL SPEED (f2 = 386 RP'I) 


106 


W (OUT-OF^-PLANE BENDING) 
V{ INPLANE BENDING) 



r/R 


(a) First Mode (Frequency X^/ft = 0.827) 

FIG. 12 MODE SHAPES FOR BOEING ROTOR AT X = 0.7 
AND ft = 386 RPH 
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FIG 



FIG 



BLADE COLLECTIVE PITCH AT 75% SPAN (9 ?5 ) 
(a) Collective Mode 


FIG. 13 BLADE NATURAL FREQUENCIES , NONROTATING, 
FOR BELL ROTOR 
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3rd (OUT-OF-PLANE) 



q 20 40 DEGREES 

BLADE COLLECTIVE PITCH AT 75% SPAN (0 <75 ) 

(b) Cyclic Mode 


FIG. 13 CONCLUDED 



NATURAL FREQUENCY PER REVOLUTION, X 



(a) First Collective Mode 

FIG. 14 BLADE ROTATING NATURAL FREQUENCIES FOR BELL ROTOR 
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NATURAL FREQUENCY PER REVOLUTION, X 



ft = 229 RPM 


ft » 350 RPM 


ft = 458 RPM (NORMAL ft) 


ft = 550 RPM 


INFLOW RATIO, v/|ft|R 
(b) Second Collective Mode 


FIG. 14 CONTINUED 
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NATURAL FREQUENCY PER REVOLUTION 



FIG. 14 CONTINUED 
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INFLOW RATIO, V/|a|R 
(a) Collective Modes 


FIG. 15 BLADE NATURAL FREQUENCIES FOR BELL ROTOR AT 
CONSTANT ROTATIONAL SPEED (0 = 458 RPM) 
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NATURAL FREQUENCY PER REVOLUTION , X.T c V|Q| 



INFLOW RATIO , V/| £2 1 R 
(b) Cyclic Modes 
FIG. 15 CONCLUDED 
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W (o) COLLECTIVE OUT-OF-PLANE BENDING 
V COLLECTIVE INPLANE BENDING 



(a) First Mode (Frequency A^°^/|(lj = 1.31) 

FIG. 16 COLLECTIVE MODE SHAPES FOR BELL ROTOR AT 
A = 0.7 AND = 458 RPM 
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4 . 93 ) 
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(d) Fourth Mode (Frequency = 10.6) 

FIG. 16 CONCLUDED 
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Second Mode (Frequency = 1.34) 

FIG. 17 CONTINUED 
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(c) Third Mode (Frequency A^ C / , '|ft| = 4.35) 


FIG. 17 CONTINUED 
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y VERTICAL BENDING 
C CHORDWISE BENDING 


<$, TORSION 



FIG. 19 MODE SHAPES FOR BOEING WING 
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d) Fourth Mode (15.9 Hz) 
FIG. 19 CONTINUED 
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VERTICAL BENDING 
CHORDWISE BENDING 
TORSION 



0.5 1.0 

WING STATION y/L 
First Mode {2.65 Hz) 


MODE SHAPES FOR BELL WING 




0 0.5 1.0 


WING STATION y/L 
(b) Second Mode (4.72 Hz) 

FIG. 20 CONTINUED 


L3G 
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IMAGINARY 



(a) Bell Proprotor 

FIG. 21 INFLUENCE OF ADDITION OF HIGHER DEGREES OF FREEDOM ON EIGENVALUES 
COMPARED WITH BASIC NINE DEGREES OF FREEDOM 
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IMAGINARY 


REAL 


O 9 DOF 
O 18 DOF 


Q* (? +D 

DO 


Q°(C) 

fid 


Q x (5-1) 

CD 


i 


- 0.2 


C>2 (8+1) 

O o 


i 2.0 


Q°(8) a 3 <p) 

o g <5* 



q 2 (e-n 

<p 

-o.i 


a 2 (q 2 ) tfi 
a x (e3 i) ® 


(b) Boeing Prop rotor 


FIG. 21 CONCLUDED 
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(a) Bell Proprotor 

FIG # 22 EIGENVALUES FOR 18 DEGREES OF FREEDOM 


i 
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(a) Bell Proprotor Eigenvalues 

FIG. 23 INFLUENCE OF BLADE MODE SHAPE TYPES ON PROPROTOR SYSTEM EIGENVALUES 
AND DAMPING RATIOS FOR COUPLED ELASTIC MODE SHAPES, COMPARED WITH 
THOSE BASED ON RIGID-BODY MODE SHAPES AND UNCOUPLED ELASTIC MODE 
SHAPES 
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(a) Collective Flapping Motion 0 Q Response to 
FIG. 25 FREQUENCY RESPONSE OF BELL ROTOR TO VERTICAL GUST u G 
INPUT AT FREQUENCY to 
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FIG. 26 FREQUENCY RESPONSE OF BELL ROTOR TO LATERAL 
GUST v_ INPUT AT FREQUENCY u 
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FIG. 26 CONTINUED 
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(a) Collective Flapping Motion e Q Response to w Q 
FIG. 27 FREQUENCY RESPONSE OF BELL ROTOR TO LONGITUDINAL GUST w G 
INPUT AT FREQUENCY (i) 
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FIG. 27 CONTINUED 
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(d) Collective Lagging Motion Response to w Q 
FIG. 27 CONTINUED 


177 



FIG. 27 CONTINUED 
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FIG. 27 CONTINUED 
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FIG. 27 CONTINUED 
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(i) Wing Torsion p Response to w G 
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(a) Collective Flapping Motion 8 Response 

FIG. 28 FREQUENCY RESPONSE OF BOEING ROTOR TO VERTICAL GUST u Q 
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FIG. 28 CONTINUED 
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FIG. 28 CONTINUED 
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FIG. 28 CONCLUDED 
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PIG. 29 CONTINUED 
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FIG. 29 CONTINUED 
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FIG. 29 CONTINUED 
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FIG. 29 CONCLUDED 
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(g) Wing Vertical Bending q 1 Response to w G 
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FIG. 32 FREQUENCY RESPONSE OF BOEING ROTOR TO LONGITUDINAL 

CYCLIC PITCH 0. INPUT AT FREQUENCY 0) 
is 


212 




APPENDIX A 


DETAILED DERIVATION OF THE EQUATIONS OF MOTION 


A. 1 Blade Motion 


The blade coordinate system is described in Section 2, and 
in Figs. 2 and 3. The position vector ^defining an arbitrary 
point on the r-axis of the nth blade with respect to the inertial 
x-y-z system becomes 


V? = vt + 


*r 

I 




l 

0 

o 


cjoU. 


COO 


coo v> p 




0 

1 

o 




caov% 


P J 




I 


uT w 


(A. 1) 


where is the position vector which expresses translation at the 
wing tip as 



(A * 2 ) 


The square of the velocity at an arbitrary point on the r- 
axis is expressed as 




= (r-nf + Cr . 1 + + (flilf + 

+ r 2 i + r x i>r + (£.*+ ** «<>*+,) 

- a r-n. r*. am. 4 » - 2&irn v’p ■*^•'-4'.. 
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+ 2 r\a Vr - Vp - 2 .rZX O p UT n A^^ K 

- mOx r x e««<l» M - 2 .^- 0 .^^ c**^ 

- 2 r x au - 2 . ft. Cr^ Cj, ^.4> n + 2 \ruu.^ 

+2ft. f“*.0 P - ar r*. O r a£~4u 

- 2r"il vV V\ <*© 4*. - 2 ft /* O r O p a**- 4**. 

- 4 lT i n vV- O r aU. «po< 4 > h 

- 2 V' 2 !! v>p O p -2 ft /■JlvvC'p "*»4* n 

-ift v-n c« 4n - 2 vA i*r K vjj cot) c|> K 

- 2 ftr 0^ Oj cw^* 4 2 ftyA vV vjj A^.4» K 

- 2 ft. 0^ ^4^ +a&sxtfUv>«j 

- 2 V\fL V\j Og A^4* n Cffo^ 4 2 Ukfi 

- 2.V uJi. Op 4- 2 V* Ol» ^ x>4n. 4« 

- 2 y 0* vf Cirqtjv + 2 rilWp fz /•*■*»- 4* 

+ 2ir r z v^ + 2^X1. vj tr z >4^.4^ 

- 2r*n V> f O p >^.4, «ro4 n 

- 2 r*OijO r AM^.^«ro 4 ' K 

+ 2^A Op ^4 * n (a. 3) 

The total kinetic energy T g for N blades is expressed as 

T E - -i- X\ ^ M ******* *** (A. 4) 

where ro is the blade mass per unit length. 
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The potential energy U E of N blades is 

Ue ° ^ J [.C 5 *)* + 

+• [(ei). + fsi) t -=*6.] (-4 t?) 

+ leer), -(«>.) (A . 5) 


The virtual work done by the centrifugal forces 6W CF is 
described as 


S ■= Z £ [-fr( T f£) ** + -y( T #) Sir "]^ <»■« 


A. 2 The Work Done by Blade Aerodynamic Forces 

Aerodynamic forces acting on the nth blade (?) n are derived 
in the z-y-z coordinate system from the and Pg components 
as 



1 

1 

*■ 

* 

-J 


" - P© " 

* 1 


£© «*®4u 
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1 

1 


Fz 


(A. 7) 


Then the virtual work done by the aerodynamic forces of all blades 

<$W is described as 
AF 


— /L» 



(A. 8) 
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A. 3 Lagrange's Equations 


In addition to the above, the kinetic energy and potential 
energy of the wing derived from the simple beam theory are append- 
ed. 

Then, Lagrange's equations are given by 

d ( ZT \ _2>T . dU_ e n. (a.9) 

it U{J su 

where T, U and are respectively, the total kinetic energy, po- 
tential energy and generalized forces of the system, and 

r 2y ^ W, (A. 10) 

Pw 

Finally, the above represents the derivation of Eqs . 2.4, 
2.5, 2.6, 2.7, 2.8, 2.10, 2.11, 2.12, 2.13, and 2.14. 
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APPENDIX B 




/ 

STRESS ANALYSIS OF THE ROTOR BLADE AND THE WING 




In this appendix, the method of deriving the stresses in 
the rotor blades and in the wing due to the gust are briefly 
stated. If the deformation of an elastic structure has been de- 
termined, it is a straightforward procedure to determine the 
stresses corresponding to this deformation. The mode acceleration 
method discussed in Ref. 15 is very suitable because the stresses 
can be determined accurately and directly when the deformation has 
been computed in terms of displacements of normal modes. 

Equation 4.57 is written here again in matrix form as 

CA1H +[8]{x| + tcH*! =LD]{el + M (b. 1 ) 

where [A], [B] and [C] are square matrices to define the coeffi- 

cients of equations derived from mass, stiffness and aerodynamic 
forces, and the [D] matrix is the excitation. The generalized 
coordinates {x} and excitation input (e) are 



{ B. 2) 


The static force {F} from the lift or drag of the wing and rotor 
in steady flight may be included if the total stresses are required. . 
When only the additional stress due to the gust or control input 
varying with time is required, the {F} matrix may be dropped. 


217 



In general, the stress at a particular point p in the rotor 
wing structure is expressed by 

KI = [VlW <B - 3 ’ 

where the [A 1 matrix is a constant matrix that represents the 
p 

stress at the point p due to a unit displacement in the normal 
mode. 

The static mode displacements are given by setting {x} = 0 
and {x} = 0 in Eq. B. 1* 

= w 1 M 4 fcr l H 

Therefore, the static stress becomes 

< b - 5 > 

When the rotor and wing are vibrating, the total displacements of 
the rotor or wing can be expressed as 

M * LC]' 1 L°1 {ej + Cel -1 {f} 

‘ B - 6 > 

Substituting Eqs.B. 5 and B.6 into Eq* B.3 7 the total stress be- 
comes 

tel -frUae - [ArJ 

- [A,][c]''[B] |i} (b.7) 

This result gives the stress in the rotor or wing at any instant 
in terms of the static stress and an additional stress due to 
the vibration. 
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